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ABSTRACT 

The radiation background above the ionization edge of He II varies strongly during and after hehum 
reionization, because the attenuation length of such photons is relatively short (< 40 Mpc) and because 
the ionizing sources (quasars) are rare. Here we construct analytic and Monte Carlo models to 
examine these fluctuations, including, for the first time, those during the reionization era itself. In 
agreement with detailed numerical simulations, our analytic model for the post-reionization Universe 
predicts order-of-magnitude fluctuations in the He II ionization rate F. Observations of the hardness 
ratio He II/H I show somewhat larger fluctuations, which may be due to more complicated radiative 
transfer effects. During reionization, the fluctuations are even stronger. In contrast to hydrogen 
reionization, our model predicts that regions with strong He II Lya forest transmission should be 
reasonably common even during the beginning stages of reionization, because of strong illumination 
from nearby bright quasars. Partly because of this, the mean ionizing background does not evolve 
strongly during and after helium reionization; it is roughly proportional to the filling fraction of He III 
regions. On the other hand, regions full of He II and also "fossil" ionized regions that contain no (or 
few) active sources appear as strong IGM absorbers. Their presence exaggerates the evolution of the 
hardness ratio, making it evolve more strongly than naively expected during the reionization era. 
Subject headings: cosmology: theory - intergalactic medium - diffuse radiation 



1. INTRODUCTION 

To the vast majority of baryonic matter in the Uni- 
verse, the most important radiation field is the meta- 
galactic ionizing background. As a result, a great deal 
of effort has gone into measuring its properties. During 
most of cosmic history, the background evolved relatively 
slowly, but there were two major exceptions: the reion- 
ization of hydrogen (at z ^ 6) and of helium (with dou- 
ble ionization occurring at z 3). During each of these 
episodes, the intergalactic medium (IGM) underwent a 
phase transition and the Universe became (mostly) trans- 
parent to the relevant ionizing photons, allowing the 
high-energy radiation field to grow rapidly. 

Recently, these events have received a great deal of at- 
tention in both the observational and t heoretical com- 
munities (see, e.g., reviews by Barkan a fc Loebl l200l 
Ciardi fc Ferrarall2005t [Fan et al.li2006i : iFurlanetto et all 



tic fluctuations in the galaxy density, and radia- 
tive transfer effects. Such variations are relatively 
small for the hydrogen-ionizing background (typically 
only a few p ercent of the r nean v alue at z ~ 3; 
ZtoI Il992allbl: iFardal &: ShuU fl99l iMeiksin fc Whitg 



2003t ICroft I2O O4') . principallv because the mean free 







ith ot tiydr ogen-ionizmg photons is extremely large 
MOO Mpc; iMadau et all 119991 : iFaucher-Giguere et all 



20061 about hvdrogen reionization), and our understand 
ing of the evolution of the ionizing radiation field has 
become increasingly sophisticated. However, most such 
studies still treat the ionizing background as spatially 
uniform. 

During reionization itself, this assumption is obviously 
wrong, because some regions are exposed to strong ion- 
izing radiation while others remain neutral with no lo- 
cal illumination. Moreover, even within the ionized re- 
gions, differences in the effective "horizon" to which ion- 
izing sources can be seen induce large spatial fluctuations 
in the background. Recent models have begun to ad- 
dress the variations expected during hydrogen reioniza- 
tion jMesinger fc Diikstral 120081: iFurlanetto fc Mesingeil 
|2009| ). demonstrating that they are crucial for interpret- 
ing measurements of that era and for understanding ra- 
diative feedback processes. 

These fluctuations will persist even in the post- 
reionization era because of source clustering, stochas- 



l2008f ) and because sources (galaxies) are relatively com 
mon. Fluctuations are larger at z > 4, when the mean 
free path is relatively small and galaxies are more highly 
clustered (e.g.. Meiksin fc White 2003 ). These studies 
concluded that such fluctuations are only significant if 
sources are rare and only measurable with future high- 
precision observations. Thus the usual assumption of a 
spatially constant background is reasonable. 

On the other hand, fluctuations in the helium-ionizing 
background have received less theoretical attention 
(though see IFarda l et al.]fl998l: iMaseUi fc Ferraral [20051 : 
iBolton et a l. 2006; iMsiksml l2007t). Nevertheless, they 
are both larger and more observationally relevant than 
those for hydrogen. One reason is that the IGM absorbs 
helium-ionizing photons more strongly than hydrogen- 
ionizing photons, leading to shorter attenuation lengths 
and larger fluctuations. Second, He II has an ioniza- 
tion potential of 54.4 eV, a sufficiently high energy that 
quasars (with hard spectra) are required to ionize it. As 
such, the sources are quite rare, and themselves quite 
variable, implying large random fluctuations in the back- 
ground. 

Moreover, helium remains singly-ionized until z ~ 3, 
a regime that is relatively easy to observe. A wealth 
of data now suggests dramatic evolution in the prop- 
erties of intergalactic helium at about this time. The 
strongest evidence comes from far-ultraviolet spectra of 
the He II Lya forest along the lines of sight to sev- 
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eral bright quasars at 2: ^ 3 (iJakobseii et all |1994|;^ 



Davidsen et al.1 Il996t lAnderson et alT 19991: Heap ct al. 
2000 HSmette et al.ll2002t iZheng et a l.ll2004bl:IShull et all 



2004 iReimers et alj" 



_ _ ■ ■ , .200 4. 2005; Fechner et all 120061: 

Fechner fc Reimersll2007( ). which show a rapid increase 
in Lya absorption by He II around that time - ahhough 
with substantial opacity fluctuations observed along sev- 
eral lines of sight ("Anderson et al."1999|: lHeap et allboOOl : 
iSmette et al. 2002; Rcimers ct al. 20Q1). These are most 
likely a direct effect of inhomogeneities in the helium- 
ionizing background during (or just after) reionization. 

Indirect measurements of the metagalactic ionizing 
background, which hardens as helium is reionized and 
the IGM becomes transparent to high-energy photons, 
suggests similar fluctuations after reionization. In par- 
ticular, measurements of the He II/H I ratio suggest 
that the helium-ionizing background fluctuates strongly 
at 2; ~ 2.6, with nearly order-of-magnitude spatial 
varia tions on scales spanning a few to a few tens of 
Mpc (ShuU ct al. 2004: Zheng ct al."2004bl: lFechner etlHI 
I2OO6; Fechner fc Rcimers _2007.). Such strong fluctu- 
ations must be due to some combination o f source 
count variations and radiative transfer effects. iMeiksinI 
(|2007f) showed with an analytic model that the broad 
distribution of quasar luminosities, together with their 
sparseness, a. c count s for much of the observed variation. 
iBolton et al.l ([lOOB) also examined this regime with nu- 
merical simulations of quasar ionizing radiation, finding 
that the former effect could account for much, but not 
all, of the observed variance. However, these simula- 
tions were limited to a relati vely small box (le s s than 
a full attenuation volume). iMaselli fc Ferraral ()2005f ) 
found (again using numerical simulations) that radiative 
transfer effects can also cause substantial variations, al- 
though they were hampered by the assumption of spa- 
tially uniform ionizing sources. ShuU et al. (2004) iden- 
tified another contributing factor: the broad distribution 
of quasar spectral indices, which directly affects the ion- 
ization rate. 

Other observations of the He II/H I ratio suggest 
rapid time evolution as well (jHeap et al.ll200dl ). which 
is qualitatively consistent with helium reionization oc- 
curring at z ~ 3: we would expect that the patchi- 
ness of the reionization process would lead to especially 
strong fluctuations between the bubbles of fully-ionized 
helium and the remaining He II lying betwe en them 
(jFurlanetto fc OhI l2008bl : iMcQuinn etall [2009) . How- 
ever, because of the difficulty of properly simulating he- 
lium reionization, these observations have not yet been 
quantitatively evaluated in light of modern reionization 
model s. The only previous study, bv iTittlev fc MeiksinI 
(|2007f ). showed that the complex radiative transfer dur- 
ing He II reionization could create wide variations in the 
hardness ratio. 

These evolving variations should also be visible in 
metal lines. For example, the ionization potentials of 
Si IV and C IV straddle that of He II, so their ratio should 
evolve during helium reionization. iSongailal {1998, 20Q5.) 
found a sudden break in that ratio at z 3 (see also 
iBoksenberg et aLll2003l) : modeling of the ionizing back- 
ground from optically thin and optically thick metal line 
systems also shows a significant hardening at z 3 
(|Vladilo et al.l[200l lAgafonova et al.ll2005L I200I . How- 
ever, other data of comparable quality show no evidence 



for rapid evolution ()Kim et al.ll2002l : lAguirre et al.ll2004f) . 

These contrasting conclusions again suggest that the ion- 
izing background may itself be fluctuating strongly be- 
tween different lines of sight. 

Here we study spatial and temporal variations in the 
helium-ionizing background using a combination of ana- 
lytic and simple Monte Carlo models. Moreover, we ex- 
amine both the relatively simp le post-reionization limit 
(as in iMeiksin fc White! l2004l ) and the behavior dur- 
ing reionization, when the fluctuations may be much 
stronger. Extending the models to the reionization epoch 
allows us to predict the distribution of the hardness ra- 
tio, and hence the observability of the He II Lya forest, 
during that era. We will show that, in contrast to con- 
ventional wisdom, substantial transmission will remain 
throughout the bulk of helium reionization. 

We describe our methods to compute the distribution 
of the amplitude of the ionizing background in Sj2] and [3l 
We describe our post-reionization results in |j4]and those 
during reionization in SjS] Finally, we conclude in Sj6l 

In our numerical calculations, we assume a cosmol- 
ogy with = 0.26, f^A = 0.74, = 0.044, = 
lOOft. km s-i Mpc-i (with h = 0.74), n = 0.95, and 
( Ts = 0-8i consistent with the most recent m easurements 
(|Dunklev et al.ll2008l: [Komatsu et al.ll2008D . Unless oth- 
erwise specified, we use comoving units for all distances. 

2. METHOD 

We wish to compute the probability distribution /(J) 
of the angle-averaged specific intensity of the radiation 
background at a given frequency (typically the ionization 
edge), J. Its overall distribution in the IGM depends on 
four basic parameters. The first two are the number 
density and luminosity distribution of ionizing sources 
(quasars in this case), parameterized by their comoving 
luminosity function <i>(L, z). The third is only relevant 
during reionization: the size of the local He III bubble, 
R. If R is finite, the IGM surrounding the bubble will 
absorb ionizing photons from sources outside the local 
region. The last parameter is the attenuation length for 
ionizing photons within each He III region, rg. We shall 
discuss To in more detail in SjSj for now we treat it as a 
constant. 

2.1. The Quasar Luminosity Function 

We use the recent estimate of the i?-band quasar lu- 
mino sity function ^ ( L, z) o ver a broad range of redshifts 
from iHopkins et al.l (|2007f) . who convert to a J5-band 
luminosity function using the observed column density 
distribution of quasars selected in the X-ray (thus ac- 
counting for obscured sources) . This luminosity function 
is consistent with earlier estimates, but it does have a 
significantly flatter faint-end slope at higher redshifts, 
d riven b y recent measurements of faint z 3 quasars 
(iHunt et al. 2004: Cri stian i et al. 2004: .Fontanot et al.1 
120071: iBongiorno et al.ll2007t iSiana et all l200/). 

At a given redshift, we define the dimensionless lu- 
minosity function (j){x) — ^{xL^,, z)/ni, where n.i is the 
total number density of ionizing sources, 

n, = J '"^^ dL^L,z), (1) 

X = L/Li,, Li, is the mean luminosity of the sample,^ 

^ Note that this is not necessarily the same as the parameter 
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and the integration ranges from minimum to maximum 
imposed luminosity cuts, Lmin and Lmax- The bright 
end of the luminosity function is sufficiently steep that 
our results are insensitive to the latter; the former is im- 
portant because a large population of faint quasars pro- 
vides a floor below which J cannot fall. We will typically 
use Lmin — 10*'^ ergs^^, slightly below the observational 
limit at z ~ 3. 

We must next convert from these _B-band luminosity 
functions to the photon energies of interest to us {E > 
54.4 eV): this requires a template for the spectral energy 
distribution of quasars. We use 

( 2500 A < A < 4400 A 

oc <^ 1050 A < A < 2500 A (2) 

A < 1050 A. 

At A > 1050 A , this template agrees with 
that of iM adau et al.l (11999 ^: other templates (e.g., 
ISchirber fc BuUockTioOl disagree in detail but do not 
affect our conclusions, given the uncertainties. Most im- 
portant for us is the far-ul t raviol et spectral index a. At 
low redshifts, iTelfer et al.l (|2002| ) find a wide variety of 
quasar spectral indices in the extreme ultraviolet, with 
a mean value of (a) ~ 1.6. This is slightly harder than 
the estimate of lZheng et al.l (|1998[ ). who foun d (a) ^ 1.8, 
but sh allower than the (a) « 0.5 estimate of lScott et al.l 
()2004[) from z <l quasars. All these studies suggest sub- 
stantial intrinsic source-to-source variance in a, which af- 
fects the luminosity function above the helium ionization 
edge. To model this, we assume that the probability dis- 
tribution of a is Gaussian in the range a G (0.5, 3.5), with 
mean a — 1.5, an d zero elsewhere; t his provides a rough 
description of the ' Telfer et al.l (I2002D data. We then con- 
volve this distribution with the ms "eTall (I200I lu- 
minosity function for our calculations. The variance in 
a slightly moderates the break near L^, and makes the 
bright end slightly shallower (from moderate luminosity 
sources with a < 1.5), but it does not have a large effect 
on our results. 

Note that our model for /(J) therefore includes self- 
consistently the variations in a among the quasar pop- 
ulation. However, when we compare to observations we 
must further transform J to the ionization rate F, which 
requires an integral over frequency and so has additional 
dependence on a: F oc (ex + 3)~ ^ in the optically thin 
hmit (see gSTl IShull et al.ll2004f ). We ignore this ad- 
ditional source of order-unity fluctuations, because our 
distributions are substantially broader anyway. 

2.2. Pure Attenuation: The Analytic Model 

We now compute /(J). We first assume that helium 
reionization has completed; this corresponds to the limit 
R 00. In that case, the attenuation length tq fully 
describes the absorption of ionizing photons. We further 
assume that Poisson fluctuations dominate variations in 
the quasar density; in that case, we can wri te the distri- 
bution analytically (|Meiksin fc White! [20031 ) 



/fl=oo(j) = 



1 



is exp 



— sNq / dx x(j){x)lmG{sx) 



sj + sNq / dx x(j){x)Ke G{sx) 



X cos 



often used in fitting the quasar luminosity function. 
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Fig. 1. — Distribution of J relative to its mean value in a fully- 
ionized IGM, {J)ionizcd- T^^^ three solid curves assume 2 = 3 
and ro = 25, 35, and 55 Mpc, from widest to narrowest. The 
other curves take ro = 35 Mpc but vary our other assumptions: 
Lmin = 10^" ergs~^ (dashed curve), fixing a = 1.5 (dotted curve), 
and assuming that all sources are quasars with number density 
equal to nq{> L*) (dot-dashed curve). 



where A^o = (47r/3)nirQ is the mean number of sources 
in an attenuation volume, j = J/Ji,, J* = L^/(47rro)^, 



G{t) 



duT^{u)e'' 



(4) 



and u — e '^Jt'^. In this limit, the mean background is 
O')ji',=oo = 3A"o, or (J) oc NqJ^, oc n^L^rQ. 

Figure [1] shows some example distributions at z = 3. 
The three solid curves take tq = 25, 35, and 55 Mpc,^ 
from widest to narrowest (or from left to right in peak 
position) . The others take rg = 35 Mpc but vary some of 
our other assumptions. The dotted curve assumes that 
all quasars have a — 1.5; it is slightly narrower than 
the fiducial model because of the reduced variance in 
the far-UV luminosity. The dashed curve assumes that 
very faint quasars exist, with Lmin — 10**^ ergs^^, three 
orders of magnitude fainter than the fiducial model. This 
increases the source density by a factor of six but has 
almost no effect on /(J), because the additional quasars 
contribute only a small fraction of the total emissivity. 
Finally, the dot-dashed curve assumes that all quasars 
have L — Li, and normalizes their number density to 
n = nq{> Li,) (so that, on average, only 0.6 quasars sit 
inside each attenuation volume). In all cases / (J) is very 
similar to our fiducial model. Evidently, the shape of the 
distribution is fixed primarily by the number density of 
bright quasars. 

Note in particular that the high-J tail is nearly in- 
variant. These points correspond to regions near bright 
quasars (their "proximity zones" ) and so depend only on 
the properties of that s ingle bright source. At s ufficiently 
high J, /(J) oc J-5/2 (|Meiksin fc White! l2003l) . 

^ These values are comparable to the mean distance that a pho- 
ton at the He II ionization edge can travel from a quasar, and they 
also comparable to the mean spacing of He II Lyman-limit systems 
at 2 ~ 3; see ^for more details. 
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Aside from the different quasar luminosity function 
and at tenuation lengths, this c alcu lation is identi cal to 
that in lMeiksin fc Whitd (|2003D and|Meiksin| ()2007D . who 
also examined fluctuations in the post-reionization back- 
ground. Our results are similar to theirs, except that we 
find the low-J tail of /(J) to approach zero instead of 
asymptote to a finite value. This has no practical impor- 
tance for the observable results. 

2.3. Finite Bubbles: Monte Carlo Approach 

During reionization, we must include the additional pa- 
rameter R, the local He III bubble's r adius (or the hori- 
zon within which sources are visible). iMeiksin fc Whitd 
([2003) describe how to obtain the full distribution /(J) 
for arbitrary R and rp (building up on the solut ion with- 
out attenuation firs t provided by IZuol Il992at see also 
iFardal fc Shulilll993h . However, we have found that the 
necessary integrals do not converge well when N, the ex- 
pected number of sources per bubble, is small. Quasars 
are sufficiently rare that this is indeed the relevant limit: 
for example, our fiducial model at z = 3 has TV « I 
when R ^ 22 Mpc - and, for L^, quasars, the expected 
number is unity only when i? ~ 40 Mpc. (This contrast_s 
sharply with the hydroge n reionization case, where N 
is always extremely large: iFurlanetto fc Mesingerll2009l .) 
Fortunately, it is precisely this rare source regime that is 
most amenable to a Mon te Carlo treatment . 

We first assume, as in IMeiksin fc Whitd (|2003D . that 
the number of quasars is Poisson distributed around its 
mean value A'^ in a region (see below for a discussion of 
deterministic clustering). For a given R, we then ran- 
domly choose the number TV of quasars in the region. 
Next, for the zth such quasar, we randomly assign a ra- 
dius r, from the central point (assuming a uniform source 
density within the ionized region) and a normalized lu- 
minosity Xi. (For simplicity, we always assume that the 
point of interest is at the center of the bubble, but we 
allow the sources to sit anywhere.) Thus the specific 
intensity from all these quasars is 

N 

^ [n/roy 

where again j = J / Ji,. The mean ionizing background is 
(j) = 3iVo(l - 6-^/'^°) (IMeiksin fc Whitd[20mh . 

There is one ambiguity in this approach: the appropri- 
ate quasar number density to use when calculating N . 
The difficulty is that quasars must sit inside of He III 
bubbles. A self-consistent reionization model would gen- 
erate the bubbles around the sources, so this requirement 
would be manifestly satisfied. However, in our Monte 
Carlo model we generate a bubble of a prescribed size 
and then randomly assign resident quasars. If we simply 
used the observed quasar number density Ui to assign 
them to the bubbles, we would find a total number den- 
sity of XYiciiini'. smaller than the observed value. 

We therefore must renormalize the input source den- 
sity to <i>(L)/xHeiii in order to preserve the true total 
quasar abundance and emissivity. This effectively de- 
mands that quasars cluster in such a way that they are 
(uniformly) overabundant in ionized regions and absent 
in the rest of the IGM. This does not accurately describe 
real quasar clustering, which is scale-dependent and sen- 
sitive to the details of the quasar host population, but it 
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Fig. 2. — Distribution of J relative to its mean value in a fully- 
ionized IGM, (J)io-nizcd- curves assume 2 = 3, ro = 35 Mpc, 
5;Hcni = 0-5i and Lniin = 10'*'' ergs~^. The thick curves take 
R = 15, 25, 35, 50, 70, and 90 Mpc from left to right in their peak 
locations (or dotted, short-dash-dotted, long-dash-dotted, short- 
dashed, long-dashed, and solid, respectively). The thin solid curve 
shows the distribution in a fully-ionized IGM. Inset: Distributions 
for discrete bubbles with no attenuation, identical sources, and 
-/V = 1, 3, and 10 (dotted, dashed, and solid curves, respectively.) 



suffices for our purposes (especially given that stochas- 
tic fluctuations appear to dominate the topology of the 
quasa r bubbles according t o recent numerical simula- 
tions; iMcQuinn et al.ir2009f) .'^ However, when compar- 
ing distributions at different xhoII, the effective source 
density (and average clustering) does differ (even if the 
comparison is done at fixed redshift). 

We then generate f{j) for a given set of (i?, tq, z, ShcIi) 
over 10^ trials. Because each such trial is time- 
consuming, we construct a grid of distributions over R, 
with the other parameters fixed. The grid spacing is 
Ai? = O.I Mpc at i? < 20 Mpc, Ai? = 0.5 Mpc up to 
i? = 35 Mpc, and Ai? = 1 Mpc to i? = 70 Mpc. Beyond 
that, we find it sufficient to use Ai? = 10 Mpc until we 
can approximate /(j) with fR=oo{j)- For intermediate 
values of i?, we select the closest bubble size in our grid; 
we have verified that the spacing is small enough that 
interpolation is not necessary. 

Figure [2] shows some example distributions for z = 3, 
ro = 35 Mpc, and Lmin = 10''^ ergs~^. The thick curves 
take R = 15, 25, 35, 50, 70, and 90 Mpc from left to right 
in their peak location. Here, for concreteness, we take 
ShcIII = 0.5. The thin solid curve shows the distribution 
in a fully-ionized IGM (from ^T^. 

Note that J f{J\R)dJ ^ 1 in some of these cases. This 
is because there is nothing in our formalism to demand 
that bubbles have at least one source; if they are empty, 
then J — 0. In practice, this is important for smaller 
bubbles (particularly R <25 Mpc). Physically, it would 
correspond to a case where a quasar appeared, ionized 

^ This ambiguity does not occur for hydrogen reionization, 
where we can self-consistently compu te the (clustered) halo 
population inside each ioni zed bubble IFurlanetto et aLl I2004al : 
IFurlanetto fc Mesingei1[2009t) . 
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Fig. 3. — Distribution of J relative to its mean value in a fully- 
ionized IGM, (J)io-nizcd- curves assume 2 = 3, ShcIII = 0.5, 
and ro = 35 Mpc. The dotted, dot-dashed, and dashed curves take 
R = 15, 25, and 50 Mpc; the solid curves show the R ^ oo limit. 
The thin and thick curves assume Lmin = 10*^ and 10'*'^ ergs~^, 
respectively. Note that the two solid curves overlap. 



the region, and then faded away - leaving a (mostly) 
ionized region with no active sources. We will revisit 
such regions in §5.41 

As one might expect, the distributions generally be- 
come narrower as R increases, because Poisson fluc- 
tuations become less important. However, note that 
the smallest fluxes are not in the smallest bubbles: 
the minimum non-zero flux in a discrete bubble is 
Lmine'-'^/^" /{AttR)^ , which actually decreases with R. 
Th e Monte Carlo distributio ns clearly converge nicely to 
the lMeiksin fc White! (|2003D hmit for large bubbles: this 
tail is filled by regions very close to a single quasar, so 
it should remain invariant as the bubble size changes.'^ 
We note that ro has relatively little effect on /(J) when 
R is finite: it is unimportant in small bubbles (because 
g-ij/ro ^ j^-j^ large bubbles simply converge to the 
i? — > oo distribution (which does not change dramati- 
cally with To; see Fig.[T]). 

In fact, in contrast to the post-reionization limit, the 
parameter that most strongly affects the distribution is 
-^min- Figure [3] compares the distributions at i? = 15, 25, 
and 50 Mpc (as well as the R ^ oo limit) for Lmin = 10'*^ 
and 10'"' ergs"^ (thin and thick curves, respectively). 
Here we have again taken XHeiii = 0.5. 

As in the i? — > oo limit, the extra sources make no 
difference for large bubbles, because they contribute only 
a small additional luminosity. However, when Lmin = 
]^q40 ergs~^ even small bubbles are likely to contain a 
source - albeit with a low luminosity - so /(J) stretches 
to significantly smaller values. It also has a larger total 
normalization. Fortunately, this uncertainty has little 

Note that, for moderately sized bubbles, the Monte Carlo pro- 
cedure overestimates the tail's amplitude; this is likely because of 
our simplified prescription that calculates the ionizing background 
at the center of each bubble; in reality, quasars can be up to two 
bubble radii away from the point of interest. This has no effect on 
our final results. 



real impact on our results, for two reasons. First, as we 
will see below the integrated hard photon background 
from quasars outside the bubble will provide a minimum, 
nearly uniform ionizing background of ~ 0.01-0.05 (J). 
Second, we are most interested in the distribution fairly 
late in reionization, when most bubbles are large. In that 
case, the variation with Lmin is modest. 

As an aside, the broad luminosity function is also cru- 
cial for another reason. The inset in FigureElshows /(J) 
in discrete bubbles with no attenuation if iV = 1, 3, and 
10 (dotted, dashed, and solid curves, respectively) a nd if 
all the sources are identical. In this restricted case, 
(|1992a|) has given an a closed-form solution for the dis- 
tribution. When TV ~ 1, /(J) has discontinuities wher- 
ever an additional source can be added; for example, if 
J < 2Lmin/(47ri?)^, only one source can fill the bubble. 
These discontinuities do not appear in our Monte Carlo 
distributions, because the wide allowed range of lumi- 
nosities smoothes over the features. They also disappear 
once N > 10, and when _R — > oo, even if all sources have 
the same luminosity, because additional sources at arbi- 
trary distances can fill in the gaps (see the dot-dashed 
curve in Fig. [T]). 

2.4. Quasar Clustering 

The models described above assume that quasars are 
randomly distributed throughout the Universe; of course, 
in reality they sit inside massive dark matter halos that 
cluster strongly. We can gauge the importance of the 
variance induced by this deterministic clustering by com- 
puting the typical standard deviation in_the number of 
quasar halos per attenuation volume, ~ ba{ro), where b 
is the average bias of quasar host halos and cr^(ro) is the 
variance in the dark matter field smoothed over a radius 
ro at the appropriate redshift. At z = 3, cr « (0.19,0.11) 
for ro = (25,45) Mpc in our cosmology. Recent z ^ 3 
surveys have found 6^8 for brjght quasars in the SDSS 
survey (|Shen et al.ll2007f ) and b ~ 5.5 for intermediate- 
lumin osity quasars selected via X-rays (jFrancke et al.l 
^200§). We wiU use the latter as more representative of 
the luminosity range over which most of the emission 
is produced. In that case, &cr(ro) ~ 0.6-1, and we ex- 
pect clustering to induce order-unity fluctuations in the 
ionizing background. This probably underestimates the 
importance of clustering, because clustering within the 
attenuation volumes cannot really be neglected. This 
may be especially important in the high- J tail, where an 
over-abundance of close quasar pairs may increase the 
amplitude of /(J). 

Fortunately, even a cursory glance at Figure [2] suffices 
to show that - even after reionization is complete - the 
Poisson fiuctuations in the number of visible quasars, to- 
gether with their wide luminosity distribution, provide 
much larger variations. We have verified that this is 
a reasonable approximation by convolving our fiducial 
post-reionization distribution with the underlying halo 
density distribution (which fixes Nq in each volume el- 
ement). This behavior differs strongly from hydrogen 
reionization, where the Poisson fluctuations are tiny and 
clustering is essential to estimating the fluctuations. 

3. THE ATTENUATION LENGTH 

One of the crucial inputs to our formalism is the at- 
tenuation length of ionizing photons, ro. Here we will 
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describe how we estimate this length scale; we will fol- 
low the method originally presented in iFurlanetto fc Ohl 
(|2008br ) and refer the reader there for more details. 

To estimate the attenuation length around an individ- 
ual quasar, we need to find the distance at which a system 
with T ~ 1 typically lies. We therefore need to associate 
overdensities with column densities in order to tr ack the 
total amount of absorption within each system. ISchavd 
(|2001f ) has shown that H I Lya forest absorbers can be 
accurately modeled by assuming their physical scale to 
be comparable to the Jeans length Lj Cg/ {GpY^^ . We 
will use the same approximation for helium, so that an 
absorber with overdensity IS. — p/ p has column density 
^Hcii « Ljx-ticiin-nc- Then 



100 



HoII 



1.8x lO^^A^/^y-o.: 
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(6) 

where we have assumed photoionization equilibrium, 
used the case- A recombination rate, and T = IQ^T^ K 
is the IGM temperature. A system will become optically 
thick at a frequency v when = (Ti/A^hcII = 1. At a 
distance R from a quasar with _B-band luminosity Lb, 
this requires an overdensity 



A,: « 64r; 



2/15 



1012 L, 







l + z 



-5/3 



R J V 4 

at the ionization edge. Higher-energy photons can 
penetrate even denser systems than implied by this 
simple model (see tj3.1l below), but the photoioniza- 
tion cross section (Ji, — cTiyHdi (^HoIi/'^)'^ (with (Jvu^^i ~ 
1.91 X 10~i^ cm^), so the instantaneous ionization rate 
(the most important quantity for us) is dominated by 
low-energy photons. However, effects that accumulate 
over time - such as heating - depend more on the 
hard photons (IAbe]_fc Haehneltlll999tlBolton et"al1l2008l : 
iMcQuinn et al.ll2009D . 

Systems with A > A^ are the "Lyman-limit systems" 
that determine rg. In the limit of a uniform ionizing 
background, ro will be the mean separation of such over- 
densities; this is the relevant c alculation during hydro- 
gen reionization, for ex ample (IFurlanetto fc Ohl l2005t 
iChoudhurv et all 120081 : IFurlanetto fc Mesingeii I2009D . 
However, we have already seen that the rarity of bright 
quasars induces large fluctuations in the helium-ionizing 
background, especially when the ionized bubbles are still 
separated by neutral walls. In this case, most points in 
the IGM are illuminated by only one (or at most a few) 
quasars. Then it is better to compute the maximum dis- 
tance from each such quasar that an ionizing photon can 
travel. 

Thus, to obtain tq, we simply find the point at which 
the mean separation of the neutral blobs, Ai(A,i), equals 
the distance from the quasar (note that \i decreases 
with distance from the quasar because its ionizing in- 
tensity decreases) ; this is the distance at which a photon 
is likely to have encountered a system neutral enough 
to absorb it. To compute the mean separation between 
IGM patches with A = A^ , we employ the IGM density 
distribution from iMiralda-Escude et all (|2000( l. as well 
as their prescription for Ai(Ai), which provides a good 
fit to numerical simulations at z 2-4. More recent sim- 
ulations show deviations from this fit at high densities. 




IQi' 10'= 
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Fig. 4. — Attenuation length as a function of i?-band quasar 
luminosity, assuming that each quasar is isolated. The solid, long- 
dashed, and short-dashed curves are for 2 = 2, 3, and 4, respec- 
tively. The filled triangles mark the luminosity- weighted mean (ro ) 
for each redshift. The dotted curve shows the maximum radius of 
a He III bubble around an isolated quasar, neglecting recombina- 
tions and assuming a lifetime of 10^ yr. 



SO our results should only be viewed as a rough guide 
(|Pawlik et al.ll2009t iBoTton fc BeckeHl2009 f) ; for this rea- 
son, we examine a relatively wide range in attenuation 
length. Moreover, we ignore a few other important ef- 
fects, like the accumulated photoelectric absorption of 
the low-density IGM; see below for a discussion of these 
problems. 

Figure |4| shows ro(L) for quasars at 2; = 2, 3, and 
4 (solid, long-dashed, and short-dashed curves, respec- 
tively). The filled triangles indicate the luminosity- 
weighted mean (rg) across the entire quasar population. 
We find that (ro) « 35-38 Mpc over this redshift range; 
the increasing dumpiness, increasing mean luminosity, 
and decreasing mean IGM density roughly cancel each 
other out. For reference, the dotted curve shows the 
maximum size of the ionized bubble surrounding an iso- 
lated source (assuming a quasar with a lifetime of 10'' yr 
shining into He II) , computed by equating the total num- 
ber of ionizing photons to the total number of He II ions 
in the bubble (thus ignoring recombinations). Even if 
quasars live for 10* yr, ro will still be comparable to the 
radius of each quasar's ionization zone. Thus, we expect 
attenuation to be relatively unimportant early in reion- 
ization but to become more significant as quasars are 
born into larger pre-ionized regions. 

Once the ionizing background becomes more uniform, 
we must include the accumulated background from more 
distant sources when estimating where neutral absorbers 
can appear. By relating the helium absorbers to the 
hydrogen Lya forest, one can show that a uniform 
background yields Ahc ~ 6.6, 12, and 30 Mpc at z = 
4, 3, and 2, with at least a factor of two uncertainty 
from the amplitude of t he helium-ionizing background 
(jFurlanetto fc Ohl[2008bl ). This suggests that the tran- 
sition from single-source to uniform is probably gradual 
and relatively smooth. Our estimates also compare well 
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to others in this regime: for example. iBolton et all (|2006f ) 
take ro = 30[(1 + z)/4]"^ Mpc. 

With the luminosity function and attenuation length in 
hand, we can now perform an important self-consistency 
check: do these parameters reproduce the observed mean 
He Il-ionizing background, often measur ed through the 
hardn ess ratio (see SJH below)? With the iHopkins et al] 
()2007D luminosity function and Rq ~ 40 Mpc at z = 2.5, 
and with (a) = 1.5, we find T « 6 x IQ-^s g-i -jhi s 
is comparable to the simulations of iBolton et al.l (120061 ) . 
which found reasonable agreement to the observations 
when a fluctuating ionizing background is included. It 
is certainly well within the errors of the measurements 
(including that of (a), which causes about a factor of 
two uncertainty in the total emissivity). We consider this 
agreement adequate, but note as well that when compar- 
ing to the observations we typically normalize the mean 
background to the observations, which is equivalent to 
adjusting (a). 

In detail, our procedure ignores a few important effects. 
First, we overestimate by up to a factor ^ 2 because 
of the accum ulated photoelectric abs orption of optically 
thin systems (jFurlanetto fc Ohll2005f l. Second, we have 
computed the mean free path at an arbitrary point in 
the IGM, whereas we are actually interested in the mean 
free path as seen by a qu asar, which most likely sits in an 
overdense region (seealsol^TL^IlOOl; \M varez fc Abell 
l2007t iLidz et al.l 120071) . This will also cause us to over- 
estimate the mean free path. However, even at z ~ 6 
(where these massive halos are much more rare) the en- 
vironments typically approach the mean density within 
< 20 Mpc of the quasar, so it should not be a large ef- 
fect. Finally, we have ignored higher energy photons, 
which can travel much farther. We will now consider 
their accumulated background. 

3.1. The Hard Photon Background 

To this point, we have assumed that, during reioniza- 
tion, all the incident ionizing radiation on a point comes 
from quasars within the local ionized volume, i.e. the 
He II in between fully-ionized bubbles absorbs all the 
radiation from external sources. However, high-energy 
photons can propagate large distances, even if all of the 
helium remains singly-ionized: for a uniform IGM, the 
mean free path is 

A„ = AndirpitriTT ( I I — —I , (8) 



, i^HcII / V 4 

where Aodgc = 0-72 Mpc is the mean free path at the ion- 
ization edge. High-energy photons will therefore create a 
diffuse, nearly uniform background that provides a lower 
limit to J, even before reionization is complete. 

To estimate this limiting value, we first assume that 
all the ionized bubbles have a fixed size and a spa- 
tial number density ribub ^ a;Heiii/(47ri?'^/3). Assuming 
(naively) that these bubbles are randomly distributed, 

1/3 

their average spacing will be A£ ~ 1.6R/x^^^^^ ~ R 
during the middle phases of reionization. For a sim- 
ple estimate, we therefore assume that all photons with 
All < A„(z/thin) ~ (A^ — R) are blocked by He II regions 
in the IGM, while others traverse a distance Atj(i^) be- 
fore being blocked. Thus fthin is the minimum frequency 
for which photons arc able to traverse one of the He II 
regions between the fully-ionized bubbles. 



We now wish to compare the ionization rate due to 
these high-energy photons, (r)^}^;^^, to the mean ioniza- 
tion rate from quasars within a discrete bubble (ignoring 
all sources outside of that bubble, and ignoring attenua- 
tion within the bubble for simplicity). This latter quan- 
tity is 



f°° dv f V 
(r)bubbio = ('^('^Hcii))H / T- - — 
Jynou "-^ V^'HoII 



h{a + 3) 



<J. (9) 
(10) 



where {Jv^^ii)^ — 3iVL^(i/Hoii)/(47ri?)^ is the mean spe- 
cific intensity inside a bubble, evaluated at the helium 
ionizing edge i^hoII and ignoring attenuation ('Zuol ll992al ). 
N = {A7T/3)R^ni is the mean number of sources con- 
tained in a bubble, and where we have assumed that all 
quasars have the mean spectral index a. 

In contrast, the ionization rate from the uniform, high- 
energy background is 



(r) 



thin - ('^(^thin))7?^=, 



dv 
hv 



t'thir 



_ Au(^'thin) 
^ (•^(^thin))^^. 



R=oo 



i^HcII 
i^thin 



(11) 

(12) 



where {J {vthin)) r^^o = 3iVoL*(j^thin)/[47rA„(z^thin)]^ is 
the radiation background evaluated at the frequency 
where the IGM becomes optically thin. The factor in 
square brackets accounts for the increased mean free path 
of high-energy photons (see eq.[8]).^ Note in equation (fT2|l 
the factor (j^Hoii/^'thin)'^, which shows the suppression of 
the ionization rate at high frequencies. 
Thus we find 
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-^u('^thin) 
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(13) 



-2-2q/3 



I+q/3 



R 



a/3 



--l-a/3 
^HcII 



(14) 



The factor with Athin can be re-expressed as w 
1)""^'^, which is close to unity. We then 



(1 65"'/^ 



find that the ratio is < 0.1 for a ~ 1.5 and R = 
15 Mpc at z = 3, so long as xhcII ^ 0.6; at larger ion- 
ized fractions, the bubb les become several times larger 
(jFurlanetto fc Ohll2008bl ). so the contribution remains at 
a few percent through most of reionization. (The ratio 
in eq. [13] becomes large very near the end of reioniza- 
tion, because this simple model ignores attenuation from 
dense clumps.) 

Thus the integrated high-energy background only pro- 
vides a few percent of the total ionization rate at any 

^ By ignoring IGM clumping, we conservatively overestimate 
the background from optically thin photons: the mean free path 
may increase more slowly than eq. JHJ in the presence of dense, 
self-shielded systems that block high-energy photons relatively e f- 
ficiently (as in the hy drogen Lya forest); see iBolton et all (gOOl); 
IMcQuinn eTall ||2009|) . 
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point during reionization. This is primarily because 
Gi, oc although the radiation intensity can be 

reasonably large, only a small fraction of the photons 
actually interact. Again, effects that accumulate over 
long time i ntervals (such as heatin g ) can still be quite 
significant (lAbel fc Haehneltl 1X9991 : iBolton et all l2008l : 
iMcQuinn et al.ll2009D . 

The high-energy photons are more important after 
reionization is over, because then Vthin ^Ucii- In that 
case, assuming Poisson-distributed absorbers, we can ap- 
proximate the attenuation length as tq (x v^] /3 = 3/2 
for a column density distribution oc N~ ^^'^, a reason- 
able approximation for the H I Lva forest (jParesce et al.l 
[1980; Zuo & Phinney 1993; Lidz et al.ll2007f l. (The scal- 
ing may approach as in eq. [8] for photons well 
above the ionization ege, but again they contribute 
only a modest amount to the total ion ization rate; see 
iBoiton et al]l2008HMcQuinn et al.ll200 a) We could then 
estimate the net fluctuations by convolving /( J|ro) with 
the spectrum-weighted ionization cross-section; we will 
take a simpler approach and estimate rg at a charac- 
teristic ionization frequency. The ionization rate from 
a logarithmic frequency interval around v is vTi, oc 
h'{L^/hh')a^ro{v) oc The effective frequency 

is therefore 



Cn'^'^'^r, _ (3 + a-/?) 



(2- 



I'Hell, 



(15) 



or ro(z/cff)/?'o(i'Hcii) = {vcs/viicuY = 1-84 for /3 = a = 
3/2. Note, however, that a full radiative transfer cal- 
culation provides a much more complex spectrum, espe- 
cially when recombination radiation is included, so these 
power-law estimat es will not be particularly accurate 
m so me regimes ()Haardt fc Madaul 119961 : iFardal et all 
fl998l) . 

This increased mean free path can make a factor ^ 2 
difference to /(J), as shown by the solid curves in Fig- 
ure [TJ We note, however, that the variance in /(J) is a 
nonlinear function of the attenuation length; in particu- 
lar, fluctuations will be strongest at the ionization edge, 
which also contributes most strongly to F. Thus this 
probably underestimates the real variation; for concrete- 
ness we will use — 35 Mpc as our fiducial value below, 
comparable to the values from Figure 2] and to estimates 
based on the Lya forest. 

4. THE IONIZING BACKGROUND AFTER HELIUM 
REIONIZATION 

As a first application of this method, we consider varia- 
tions in the ionizing background after helium reionization 
is complete. In this regime, our analytic model from §2.21 
describes the distribution of /(J) and hence /(F) (once 
an appropriate average attenuation length is chosen).^ 
For example. Figure [T] showed some post-reionization 
distributions at z = 3. Regardless of our assumptions 
about the minimum quasar luminosity and the attenua- 
tion length, /(F) has a large variance. We parameterize 

^ There is one additional complication: scatter in a directly 
affe cts the ionization rate through the high-energy photons (eq. 
I15| l. over an d ab ove the additional luminosity dispersion that we 
included in ^TTI ( ShuU et al. 2004). We do not include this ad- 
ditional scatter, although it only introduces extra fluctuations of 
order unity. 



the variation by the range AF for which /(F) > /((F})/2. 
For the curves in Figure [T] with — (25, 35, 55) Mpc, 
we find AF/ (F) = (1.3, 1.15, 1.07). Order unity fluctu- 
ations are typical, and a non-negligible fraction of points 
can have much larger ionizing backgrounds. (Note, how- 
ever, that the median F is smaller than the mean.) 

This distribution should be observable via the Lya for- 
est; in particular, by combining information about the 
hydrogen and helium forests we can measure the frac- 
tional fluctuations in the ionizing backgrounds for H I 
and He II. The former is most strongly affected by pho- 
tons near the hydrogen-ionization ed ge; at z ^ 2-3 their 
mean free path is hundreds of Mpc (jMadau et al.lll999l : 
iFaucher-Giguere et al.l 120081 ) . and each attenuation vol- 
ume contains a huge number of quasars and star-forming 
galaxies. We therefore assume th at the hydrogen- 
ioniz i ng backgro und is truly uniform (|Meiksin fc Whit^ 
l200l[Croftl[200l . In that case, the He II/H I ratio serves 
as a proxy for /( F) for helium. 

The equivalent iGunn fc PetersonI ([1965) IGM optical 
depth due to He II in a parcel of gas with relative over- 
density A — pIp'is. 

3/2 



THcii = 3.55 X lO-^XHciiA 



(16) 



In ionization equilibrium with a helium-ionizing back- 



ground of amplitude F = 10 



1 



^"ri4 s- 
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this is 



aHcii 



aA(1.5 X 104 K) 



(17) 

where we have used case-A recombination under the as- 
sumption that most of the recombination photons will 
escape the local volu me and ionize denser regions (e.g., 
lMiralda-Esc^[2003h . 

Variations in THeii can therefore be traced to den- 
sity structure in the IGM, temperature variations, or 
the ionizing background. We will ignore temperature 
fluctuations for simplicity, though they may also be im- 
portant if Lya forest lines are ther mally broadened by 
high post-reionization temperatures (iGleser et al.l l2005l : 
iFurlanetto k Ohll2008al: IMcQuinn et al.ll2009f) . Densitv 
variations can be calibrated by comparison to the H I 
Lya forest, so we can isolate the effect of F by measuring 
77 ~ -/VhcIi/A^hi, which is (nearly) directly observable by 
comparing lines in the hydrogen and helium Lya forests. 
Assuming the gas is highly-ionized and in photoioniza- 
tion equilibrium, we have 

^ _ "HcIII QHcII Fhi _ ??o (F) 
nmi am F F 
We normalize the distribution of 77 by fixing its 
mean to a specified value. Unfortunately, this choice 
is not straig htforward. Theoretical estimates spa n 
in) ^ 40-80 (iHaardt fc Madaij [T99^ iFardal et al.l[T99l . 



IBolton et al.l (|2006D found that setting (r?) w 60 re- 
produces the observed optical depth in the helium 
Lya forest when a fluctuating radiation field is in- 
cluded. Observations show a wide scatter in the val- 
ues of T] in individual syste ms, but the avera ge is ^ 45- 
110 p iuU et al. 2004; Zhen g et al.l l2004b: Fcc hner et all 
l2006n^echner fc ReimersI I2007| ). with most of the un- 
certainty due to systematics in interpreting the spec- 
tra and in how to treat outliers. For concreteness, we 
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Fig. 5. — Differential and cumulative probability distributions 
of the hardness ratio, rj, assuming that helium reionization has 
ended (top and bottom panels, respectively). The two solid curves 
take 2 = 3 and 2.5, assuming ro = 35 Mpc. The short-dashed 
curve assumes ro = 55 Mpc, respectively. The long-dashed and 
dotted curves show the "clean" subsample and full data set of 
IFechner et al.l I I2006I) . respectively. In the bottom pan el, the long 
dash-d otted curve shows the observational data from IShuII et al.l 
(120041 ) . and the short dash-dotted curve shows the simulated dis- 
tribution of lBolton et al., (i2006i) . 



will set (?7) — 80, f or easy comparison to the data of 
IFechner et al.l (|2006D . Note, however that this i s some - 
what larger than estimates from, e.g.. lShull et al.l(|2004D . 

Figure |5| shows the resulting differential and cumu- 
lative distributions of the hardness ratio in the post- 
reionization limit. The solid curves assume vq = 35 Mpc, 
at = 3 and 2.5. The variations over this redshift range 
are extremely small, because the quasar luminosity func- 
tion is not evolving rapidly {Nq ~ 13-19 from z = 3-2, 
if ro is constant). Thus any observed evolution in /(jy) 
over this interval would indicate evolution in the mean 
free path. 

The short-dashed curve shows f{r]) for ro = 55 Mpc 
at z = 2.5. Just like /(F), the hardness parameter has 
a larger variance if the attenuation length is smaller, be- 
cause Poisson fluctuations in the quasar counts become 
more important. 

The short dash-dotted curve in the bottom panel shows 
the cumulative distribution from the lBolton et al.l ()2006D 
simulations, which included discrete, clustered sources 
within each attenuation volume but did not model radia- 
tive transfer through a clumpy IGM. Their distribution 
is somewhat broader than ours (note that we have renor- 
malized the values to our choice of (77)). One reason is 
their treatment of attenuation: their simulation box was 
~ 30/1""'^ Mpc across, comparable to their assumed ro. 
To simplify the effects of absorption, they only allowed 
sources within one attenuation length to illuminate any 
given point, but they ignored absorption within that re- 
gion. This likely accounts for at least part of the dis- 
crepancy, because the added sources at large distances 
will help to damp out smaller scale fluctuations (for ex- 
ample, compare the dot-dashed and thin solid curves in 
Fig. [2]). The rest is likely due to the much more realistic 



treatment of Lya forest features in the simulations. 

4.1. Comparison to Observations 

The long-dashed and dotted curves sho w the observed 
T] distributions of IFechner et al.l (|2006f ). The dotted 
curve is their complete sample; the long-dashed curve 
is their more reliable sample of moderate optical depth 
lines, where the comparison of H I and He II columns 
can be made most reliably (these systems have 0.01 < 
Tm ^ 0.1). The agreement between our model and this 
restricted sample is quite impressive: the only signifi- 
cant discrepancy is a longer tail toward high 77 in the 
data. This may arise partly from an underestimate of 
(77) ; iMeiksin (2007), who performed a similar comparison 
(and obtained results similar to ours), found that such a 
modification would improve agreement at that end. 

On the other hand, the full data set shows relatively 
poor agreement, with significant tails toward both large 
and small rj. Other data show similar trends; for ex- 
ample, the long dash-dotted curve in the bottom panel 
show s the observed cum ulative probability distribution 
from iShull et al.l ([20?)1 ). using FUSE measurements of 
the He II forest (smoothed on 0.2 A scales) and high- 
quality optical data of the H I Lya forest along the 
line of sight to HE 2347-4342; the data span the range 
z = 2.3-2.9. They claim to measure rj in the range 
0.1-460 reliably; they have censored points within la 
of either zero or complete transmission, which does com- 
plicate the comparison.^ Zheng et al. (2004b) indepen- 
dently analyzed the same data with different techniques 
(line-fitting, as opposed to pixel methods) and obtained 
a similar distribution. 

Both of these full distributions are considerably 
broader than the theoretical predictions, f or any r cason - 
able attenuation length, and the IFechner et al.l (j2006l ) 
subsample. This is not surprising, because a variety 
of effects other than shot noise in the quasar counts 
contribute to the observed distribution, including ob- 
servational errors, clustering, temperature fluctuations, 
line broadening and peculiar velocities in the forest, ex- 
tra fluctuations from variations in quasar spectral in- 
dices, and radiative transfer effects throu gh the IGM 
(|Abel fc Haehneltl 11999': 'Bolton et al.l l2008f l. Moreover, 
especially with thefShuU et al. (2004) data, the observa- 
tions span a relatively large redshift interval, over which 
(r/) evolves significantly. 

Overall, the agreement with the "clean" data s eems re- 
markably good, as it was in the simpler model of lMeiksinI 
([2007). Shadowing and similar radiative transfer effects 
may be responsible for the high-Ty tail, which co rresponds 
to a s mall He II ionizing background (|Tittlev fc Meiksinl 
|2007| ). The good agreement with this restricted set sug- 
gests that the largest problem with the comparison to 
the full data sets arises from observational errors. 

Two other asp ects of t he da ta merit some discus- 
sion. First. iShuU et all (|2004l ) and, with less confi- 
dence, IFechner et al.l (|2006| ) detected an anti-correlation 
between A'hi and 77, so that void-like regions in the 

^ In particular, the mean of this censored data set is (rj) = 
45, smaller than our distribution. Thus the censored points lie 
primarily at very large rj, extending the observed distribution to 
the right of the plot. The differences with our models described 
in the text are qualitatively unchanged if we rescale our models to 
match this mean. 
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IGM tended t o hav e softe r ambient radiation fields . 
iFechner et aH (|2006l ) and iFechner fc ReimersI (|2007l ) 
showed that this is at least partly due to noise in 
low-column density systems and saturation in higher- 
column density systems, and they also argued that the 
anti-correlation would disappear if higher-column density 
lines were thermally, rather than turbulently, broadened. 
In any case, our model cannot address such a trend, as we 
do not include the density field in our calculations. The 
best we can say is that dense absorbers are most likely 
to be near the massive quasar hosts, which would prob- 
ably cause a weak association with low-ry lines. Naively, 
voids should be preferentially illuminated by a softer, 
more highly-attenuated radiation field because they are 
farther from the sources, at leas t if the hyd rogen-ionizing 
background itself is uniform f,Bolton et al...2006 ') . Radia- 
tive transfer, on the other hand, may cause a corr elation 
in the op posite direct i on (iMaselli fc Ferrarall2005[ ). 

Second. IShuU et all ()2004D detected substantial rj fluc- 
tuations on rather small scales (~ 2 Mpc) toward HE 
2347-4342. The natural fluctuation scale in our model is 
the at tenuation length, tq. However. IFechner fc ReimersI 
()2007f l showed that the highest signal-to-noise por- 
tions of that spectrum, and the entire line of sight to 
HS 1700-f6416 (also high signal-to- noise) , show much 
smoother variations, with ^ 33% of the IGM varying 
on scales < 6 Mpc and > 50% varying only on scales 
> 14 Mpc. The latter is not so far from our expected 
attenuation lengths, and additional effects such as radia- 
tive transfer, density structure in the IGM, and aliasing 
may explain the differences. 

5. THE IONIZING BACKGROUND DURING HELIUM 
REIONIZATION 

5.1. The Bubble Size Distribution 

We now turn to the helium ionizing background before 
overlap, when portions of the IGM are still filled by He II. 
In this case, our model requires one additional ingredient: 
the size distribution nb{R) of discrete ionized bubbles at 
a given time during reionization. 

F 'or this purpose, w e use t he excursion-set approach 
of iFurlanetto fc Oh! (|2008H) . which is based on a 
model originally developed for hydrogen reionization in 
IFurlanetto et al.l ([2004b). In brief, the model generates 
He III regions by comparing the number of ionizing pho- 
tons generated inside the region to the number of helium 
atoms; large-scale overdensities host more massive ha- 
los (and hence more quasars) and so ionize themselves 
earlier. The model uses the excursion set formalism to 
compute the scales at which any given IGM parcel enters 
an ionized region, and this is interpreted as nf,(i?). 

The two inputs for this model are xhgIII and a pre- 
scription for the helium-ionizing fluence of dark matter 
halos (essentially, a model for placing quasars in their 
host galaxies); the resulting ni,{R\x}icUi) is nearly inde- 
pendent of redshift. As a fiducial model, we will as- 
sume that halos with m > 5 x 10^^ Mq host quasars, 
and that the total fluence is proportional to the halo 
mass.® This provides a reasonable match to clustering 
measurements of luminous, high-z quasars and to the 
measured redshift dependence of the total quasar emis- 

* Note that this assumption does not affect the quasar number 
density, which we fix to the observed value. 



sivity (IFurlanetto fc Ohll2008lJ) . Physically, it assumes 
that every massive halo contains a supermassive black 
hole with a fixed radiative efficiency. 

Although the excursion-set model, which is driven 
by deterministic halo clustering, works very well for 
hydrogen rei onization ( i n comparison to numerical 
simulations: IZahn et all 120071 : iMcQuinn et all 120071 : 
iMesinger fc Furlanettoll2007l ). it is not yet clear how suc- 
cessful it will be with helium reionization, due to several 
complications during the latter era. First, quasars are 
much rarer than the galaxies that (probably) dominated 
hydrogen reionization, and the same Poisson fluctuations 
described in our model for f{J) will also dominate the 
topology of ionized gas during the early (and perhaps 
middle) stages of reionization. We mimic this by setting a 
minimum bubble size i?min during the early stages, when 
the rare quasars are sca ttered about nearly randomly 
(IFurlanetto fc Ohll2008lJ) . We set i?min = 15 Mpc, ap- 
propriate for an quasar shining for 10^ yr (see the 
dotted curve in Fig. |4|). 

Second, the relatively hard spectra of quasars allow 
some high-energy photons to propagate large distances 
through the IGM before being absorbed (as in 
These will gradually ionize and heat regions far from 
quasars. Thus it is not formally correct to divide the 
IGM into pure He II and pure He III, as the excursion- 
set model implicitly does. However, most ionizing pho- 
tons are absorbed within a few Mpc of their host bubble 
(e.g.. lMcQuinn et al.ll2009( l: so long as the characteristic 
bubble size is larger than this, our approximation is rea- 
sonable. We account for the remaining photons via the 
hard, weak ionizing background described in i )3.1l 

Third, the simple "photon-counting" model we use 
here does not properly incorporate recombinations, 
which slow the growth of large bubbles. These can be in- 
cluded in_the excursion-set formalism (Furlanetto fc OhI 
120051 l2008b). but in that case the method returns the 
distribution of the minimum of (i?, ro). For simplicity, 
we take the more transparent route of computing the 
sizes of discrete regions and then imposing attenuation 
when calculating the ionizing backgro und. 

Despite these shortcomings, the IFurlanetto fc OhI 
(2 008b) model is still the only concrete calculation of 
bubble sizes, so we will use it here with the caveat 
that the distributions are probably not quantitatively 
accurate. The crucial aspect is that it has the correct 
qualitative behavior, in which discrete bubbles gradually 
grow larger and larger, until they surpass the attenua- 
tion length (beyond which the total size matters less and 
less). Moreover, note that the bubble size distribution at 
early times is not essential in the limit of rare sources; 
for example, if each bubble hosts only one source, the 
ionizing background at interior points is independent of 
R - which only determines the minimum allowed J. 

5.2. Results 

Once ni,{R) is specified, we compute the overall distri- 
bution of the amplitude of the ionizing background via 

f{J\ro)^ J dRnbiR)fiJ\R,ro), (19) 

where we use the Monte Carlo model for R < 90 Mpc 
and equation ([3]) for larger bubbles. 
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Fig. 6. — Distribution of J during reionization, relative to its mean value in a fully-ionized IGM, (</)ionize(j- The solid curves in both 
panels show /(J) in our fiducial model, with 2 = 3, -Lmin = 10**^ erg s~^ and rg = 35 Mpc, at SHelll = 0-3, 0.5, 0.75, and 0.9, from bottom 
to top at the peak. The dot-dashed curves show the post-reionization distribution for the same parameters. Left panel: The dotted curves 
take -Lmin = 10*'' erg s~^, for the same four points during reionization. Right panel: The dotted curves take ro = 25 Mpc. 



Figure [B] shows the resulting distribution in several spe- 
cific cases. The solid curves in both panels show our fidu- 
cial model, in which z = 3, Lmin = 10**^ erg s^^ and vq ~ 
35 Mpc. From bottom to top (at the peak), the curves 
assume 5;Heiii = 0.3, 0.5, 0.75, and 0.9. For comparison, 
the dot-dashed curve shows the post-reionization distri- 
bution for the same parameters. 

In general, the tre nds are similar to those duri ng hy- 
drogen reionization (jFurlanetto fc Mesingeill2009( l. Ini- 
tially, the ionized regions are all relatively small, so any 
individual point is illuminated by only a few sources and 
the amplitude of the ionizing background is usually rel- 
atively small. As bubbles grow and encompass more 
sources, the amplitudes increase. The rate of increase 
slows once attenuation becomes important (when the 
mean bubble size exceeds rg), even t hough R actually 
accelerates its growth in this regime (jFurlanetto fc OhI 
[2008b). The distributions then match smoothly onto the 
post-reionization /(J). The final stage is the incorpo- 
ration of the remaining small bubbles into large ionized 
regions, destroying the low-J tail as these previously- 
isolated points now receive weak illumination from dis- 
tant quasars. 

However, there are important differences with hydro- 
gen reionization. First, there is less evolution at large 
J here. This is partly because the rarity of sources 
makes even the post-reionization distribution broader 
and partly because the high-J tail is created by points 
close to a single quasar, so its amplitude depends only on 
their abundance - which is fixed between all these curves. 
Thus, the transition from pre- to post-reionization is even 
smoother than for hydrogen reionization. 

Another feature that differs from hydrogen reioniza- 
tion is the normalization of /(J). If sources are com- 
mon enough that every bubble contains at least one, 
then the integral of /(J) must equal ^Heiii- How- 
ever, during helium reionization quasars are sufficiently 



rare that many smaller bubbles contain no sources. 
For example, the mean number of sources per atten- 
uation volume is only 4.2 after reionization is com- 
plete (if ro = 35 Mpc at z = 3); thus any individual 
15 Mpc bubble has only an ~ 8% chance to contain 
a quasar (although, during reionization, this increases 
because quasars can only live inside ionized bubbles). 
Thus, when ShcIII — (0.3, 0.5, 0.75, 0.9), the actual frac- 
tion of space illuminated by ionizing radiation is only 
(0.09, 0.22, 0.55, 0.81). This has interesting observable 
implications; see §5.31 

Moreover, we also find that the mean J within ionized 
regions is nearly constant across all these models: it in- 
crease by only ~ 20% from xhoIII — 0.3-1. Thus the 
average J, including both He III regions and the He II 
walls between them, is roughly proportional to XHeiii- 
This is not too surprising; we have held the comoving 
emissivity constant between these models, so J can only 
increase by allowing photons to propagate farther. We 
have rp / i?min ~ 2.3, so the mean background can only in- 
crease by roughly that factor. Provided our estimates for 
the attenuation length are reasonable, helium reioniza- 
tion will therefore not be accompanied by an enormous 
increase in the ionizing background. 

The two panels of Figure [6] also show how these results 
depend on some of the model's input parameters. In the 
left panel, the dotted curves assume Lmin — 10 erg s 
In this case, ionizing sources are much more common, 
and > 80% of of ionized regions contain active sources 
even when irneiii = 0.3. However, these additional 
sources are all extremely faint. In reality, the uniform 
hard photon background ( §3.ip will provide a lower limit 
to the ambient ionizing background of ~ 0.05 (J); Fig- 
ure [5] then shows that extremely faint active galaxies will 
not affect the observable distribution. 

The right panel of Figure [6] shows the distributions if 
ro — 25 Mpc. Perhaps surprisingly, this makes very little 
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Fig. 7. — Distribution of J during reionization, relative to its 
mean value in a fully-ionized IGM, (J) iodized- "^^^ solid curves 
show /(J) in our fiducial model, with z = 3, Z/niin = 10*^ erg s~^, 
and ro = 35 Mpc, at xhoIII = 0.3, 0.5, 0.75, and 0.9, respectively. 
The dot-dashed curves show the post-reionization distribution for 
the same parameters. The dotted curves show the same four phases 
during reionization, but allowing quasars to sit inside smaller halos 
(see text). 



Fig. 8. — Distribution of the hardness ratio, 77, during helium 
reionization. The solid curves show our fiducial model, with z = 3, 
Lniin = lO**^ erg s~^, and ro = 35 Mpc, at ShoIII = 0.3, 0.5, 0.75, 
and 0.9, from bottom to top. The dotted curves are same, except 
they allow quasars to live in smaller dark matter halos; note that 
the xhcIII = 0.3 curves overlap in both panels. The dot-dashed 
curve shows the post-reionization distribution, for the same pa- 
rameters. 



difference to /(J). One reason is that we have scaled to 
the mean value after reionization, which is proportional 
to To; for a fixed quasar emissivity, (>/)ionized 
crease with the attenuation length. However, the shapes 
of the distributions are nearly invariant, which can be 
understood by referring to Figure [1] and by noting that 
attenuation has only a minimal effect on small bubbles. 

Figure [7] shows how /(J) depends on nb{R). The solid 
curves again show our fiducial model. The dotted curves 
are identical, except they assume that quasars sit inside 
halos with m > 5 x 10^° Mq, an order of magnitude 
smaller than the fiducial model. ^ These halos are less 
strongly clustered, which shifts nb{R) toward somewhat 
smaller radii. The net effect is that more of the Universe 
sits inside of small bubbles, so /(J) shifts leftward (and 
its normalization decreases, because bubbles are more 
likely to be empty) . Evidently the details of the bubble 
size distribution do not affect the qualitative trends of 
the model. 

We should note that we have assumed a fixed red- 
shift {z = 3) and attenuation length in these figures. 
Of course, reionization will actually occur over a fi- 
nite time inte rval (although a small one, if current es- 
timates hold; iFurlanetto fc OhI l2008bl: iMcQuinn erall 
over which the luminosity function and the at- 
tenuation length itself may evolve. Figure [1] shows that 
evolution in the luminosity function has only a min- 
imal effect. Moreover, Figure [4] shows that attenua- 
tion around individual quasars - the limit most rele- 
vant during reionization - does not evolve either. On the 
other hand, once a more uniform background builds up 
(at xhoIII ^ 0.8, according to Fig. [6|), estimates based 

^ This does not affect the number density of active quasars - 
which need not sit inside every such halo. 



on the H I Lyg forest predict rp oc (1 + z) ^ or so 
(jMadau et al.lH999t [Faucher-Giguere et al.ll2008D . This 
may accelerate the evolution during the final phases of 
reionization, although the short available time interval 
makes the total change small - so we still securely ex- 
pect smooth evolution of the ionizing background onto 
its post-reionization form, ra ther than a sharp jump (c.f. 
IFurlanetto fc Mesingeil l2009 for H I reionization). 

5.3. 77 During Reionization 

The hardness parameter 77 provides one way to con- 
strain this smooth evolution. As with our post- 
reionization distributions, we can easily transform /(J) 
into f{r]) during reionization, assuming that all points 
are in ionization equilibrium (see below). Figure [8] 
shows the results. The upper and lower panels show 
the probability distribution and the cumulative distri- 
bution function, respectively. In each panel, the dot- 
dashed curve shows our fiducial post-reionization model, 
with z = 3, Lmin — 10'*'^ ergs"-^, and ro = 35 Mpc. 
The solid curves take the same parameters but assume 
ShoIII = 0.3, 0.5, 0.75, and 0.9, from bottom to top. 
Again, we have normalized these so that, in the post- 
reionization Universe, (77) = 60. 

The distributions during reionization are markedly dif- 
ferent from those after the process has completed, princi- 
pally because a large fraction of the points have ry ^ 100 
- and indeed those either entirely outside of ionized bub- 
bles or inside of an empty one have an essentially infinite 
value (see ^5A\ for a discussion of the latter) . Interest- 
ingly, the fraction of such pixels is significantly larger 
than aiHcii, because of the many points inside empty 
(but mostly ionized) regions. This exaggerates the "ap- 
parent" He II fraction that one might naively associate 
with the fraction of pixels with 77 ^ 100 and appears 
to make helium reionization a sharply defined, dramatic 
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event with this observable (even though, as we have em- 
phasized above, /(J) does match smoothly onto the post- 
reionization form, and in a regular fashion). 

This smoothness does manifest itself in /(r/) as the 
relative stability of the low-rj tail. Our models predict 
that, at least during the latter half of reionization, a 
substantial fraction of the Universe will have 77 < 100: in 
other words, regions of high transmission will continue to 
appear well into the helium reionization era. This is a key 
difference from hydrogen reionization, where the overall 
saturation of H I absorption makes studies of the Lya 
forest almost impossible before the hydrogen becomes 
fully ionized. Note that these regions do not need to 
be extraordinarily close to bright quasars; they simply 
need the local background to be comparable to the post- 
reionization value, or to have a few L^, quasars within a 
couple of attenuation lengths. 

Thus the differing character of helium reionization - 
and especially the more rapid buildup of large ionized 
regions and the dominance of rare bright sources - sug- 
gests that spectral studies of that reionization era will 
continue to provide a great deal of information about 
the IGM; we do not by any means expect complete sat- 
uration in the forest at 2; > 3. FoUowup of the many po- 
tential He II Lya forest lines of sight recently discovered 
through the Sloan Digital Sky Surve y and the GALEX 
mission should be extrem ely fruitful (iZheng et al.ll2004aL 
[2008t iSvphers et al.l[200l . 

Recall tha t the "reliable" r/ distribution from 
iFechner et al.l (|2006l ) showed a longer tail toward high 
77 than permitted in our post-reionization distributions. 
This could, in principle, indicate that this line of sight 
is probing the patchy phase of He II reionization, where 
small ionized bubbles give much bigger values of rj. How- 
ever, we would expect a very long tail if this were the 
case, rather than simply extending slightly farther than 
the post-reionization models. We therefore regard this 
interpretation as unlikely at this point, although if it 
were the case then this would indicate that He II reion- 
ization did not complete until z < 2.7. 

The dotted curves in each panel show the distribu- 
tions for a model in which quasars are allowed to inhabit 
smaller dark matter halos (as in Fig.[7|), so that the ion- 
ized bubbles are also smaller. (Note that the XHeiii — 0.3 
curve overlap in both panels.) In this case, they are more 
likely to lack active sources, so the curves lie below the 
fiducial model. In principle, this can be used to measure 
(crudely) the characteristic bubble size at a given stage 
of reionization, although at the moment it is probably 
more accurate to view the difference as the uncertainty 
in simple models. Detailed interpretations of the data 
will require careful compariso n t o simulations, s uch a s 
those of .McQuinn et all (|2009D or lPaschos et all (|2007D . 
This is particularly important toward the end of reion- 
ization, when the growing hard photon background may 
help to eliminate the last vestiges of He II, making the 
evolution during these last phases smoother than in our 
simple models. 

One possible complication is the damping wing of He II out- 
side of the ionized region, which is extremely important during H I 
reionization (Mir alda-Escudc 1998). However, the optical depth 
of He II is two orders of magnitude smaller, so the wing is much 
weaker during this time and will not interfere with most of the 
large ionized bubbles that we expect. 



Simulations are also necessary to quantify the com- 
plex effects of radiative transfer during He II reioniza- 
tion, which may br oaden the distribution substantially 
(jTittlev fc: Meiksinl 2007), including especially a wider 
tail toward high values of rj. We have also ignored varia- 
tions in a, which introduce additional fluctuations in the 
frequency-integrated F. 

We have neglected the nearly uniform, optically thin 
background in Figure [g but it provides J ~ 0.05 (J) 
everywhere. At this level, 77 ~ 1200, off the right edge 
of our plot. Measuring the growth of this weak back- 
ground will probably be extremely difficult. If the Lya 
forest lines are broadened only by turbulence or peculiar 
motions, the hardness ratio can be related to the opti- 
cal depths of pixels in t he H I and He II forests through 
r] — 4thcIi/thi (Mir alda-EscuddflQQS? ). In practice, reli- 
able measurements require 0.01 ^ t < 1, which limits the 
total measurable range of the hardness ratio to 77 ~ 1- 
400. Thus, the high-energy background is probably not 
measurable with line ratios. 

5.4. Fossil Bubbles 

We now re- visit the properties of ionized bubbles with- 
out any active sources. Above we assumed that they have 
V ^ iv) ^-nd were essentially invisible ~ even though they 
still may be highly ionized overall. These regions initially 
sat in photoionization equilibrium with their original ion- 
izing source. Once that source shut off, they fall out of 
equilibrium (in the sense of being over-ionized) and be- 
come "fossil" bubbles that recombine rapidly ~ although 
also non- uniformly - thanks to the dumpin ess of the in- 
tergalactic medium (jFurlanetto et al.ll2008f l. The ratio 
of the helium recombination time to the Hubble time is 

0-4 1^ 4 y/^^ a^(2xl04K) j^ ^^q) 



H-\z) A VI , . 
where we have assumed matter domination and a slightly 
higher temperat ure than before (because the region was 
recen tly ionized; iFurlanetto &: Ohll2008al : iMcQuinn et al.1 
[2009I) . Thus these regions will remain relatively highly 
ionized for of order the expansion timescale - particu- 
larly in the low-density voids which provide most of the 
transmission in the He II Lya forest. 

However, strong absorption in He II requires only a 
small fraction of the material to have recombined (see 
eq. 1161) . How rapidly will a typical feature become sat- 
urated in absorption? For a rough estimate, we let 
a^Hcii ~ a(r)?7.e(A)tfossii (valid when xhcII ^ 1) and find 

Thus, all regions with A > 0.1 will already be saturated 
He II absorbers just 50 million years after the quasar 
shuts off. In the simple Lya forest model of ISchavd 
()2001h . this includes all Lya forest systems with H I col- 
umn densities > 10^^ cm~'^ (or thi ^ 0.01) at z 3, 
and in simulations su ch systems cover > 9 5% of the vol- 
ume of the Universe (jMiralda-Escude et al... 2000.) . Thus 
the vast majority of the observable IGM will recombine 
quickly after their illuminating source expires. 

The recombining regions will provide a smoother tran- 
sition to high rj values, but they will typically be very dif- 
ficult to observe. Only rare voids can still provide trans- 
mission over long timescales. (Of course, all these regions 
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will eventually equilibrate with the uniform, high-energy 
background of §3.11 But, as noted previously, this back- 
ground is not strong enough to keep absorbers transpar- 
ent.) 

6. DISCUSSION 

We have examined the evolution of the high-energy 
ionizing background during and after helium reioniza- 
tion. Through a combination of analytic calculations 
and a Monte Carlo model, we have computed the prob- 
ability distribution of J, the amplitude of the ionizing 
background at the ionization edge (which is a reasonable 
proxy for the ionization rate, F). After reionization is 
over, our model requires only the attenuation length of 
helium-ionizing photons and the quasar luminosity func- 
tion (in the relevant energy range) as inputs. We find 
that bright quasars remain rare enough at z ~ 2-3 that 
/(J) is quite broad, with a fuU-width -at-half-maxim um 
comparable to (J), in agreement with iMeiksinI (|2007( ). 

Currently, the distribution of the hardness ratio 77 = 
Nucii/Nm is the best way t o mea s ure such fluctuations 
tZhcng ct al. 2004b; ShuU _^al.l l2004t iFechner et all 
I^OOfi : .Fechner fc Reimers, 20071) . The best-observed dis- 
tribution iFechner et al.l ()2006l ) is slightly broader than 
our models predict, especially toward high 77, but the 
overall agreement is good (see also Mciksin 2007). Some 
of this additional broadening is undoubtedly due to more 
complex radiative transfer than we include here as well as 
variations in the spectral indices of the ionizing sources. 
Note that there also appear to be observational bi- 
ases affecting many of the estimates in the literature. 
For example, the apparent optical depth method be- 
comes unreliable iii satur ated lines (jFechner et al.ll2006l : 
IFechner fc Reimersll2007t ). A more careful comparison of 
detailed theoretical models and the data is still needed 
in order to quantify the tension. 

During reionization, quasars sit inside of discrete ion- 
ized bubbles; the "walls" between them are full of 
He II that blocks low-energy photons from more distant 
sources. Thus f{J) depends on the size distribution of 
these bubbles. With a simple analytic model of this dis- 
tribution, we have used a Monte Carlo method (that also 
includes attenuation) to generate f{J) in this regime. 
Initially, f{J) is quite broad because of the range of bub- 
ble sizes (and because most bubbles are sufficiently small 
that a faint quasar can provide all of the local ionizing 
radiation field). As the bubbles grow larger, they en- 
compass more and more sources and the variance slowly 
declines. Eventually, the distribution matches smoothly 
onto its post-reionization form: after XHeiii ~ 0.9, the 
only significant change is the disappearance of the low- 
J tail from isolated regions. Thus we do not expect a 
sharp feature in the overall amplitude of the ionizing 
background at the completion of reionization. 

We found that the mean ionizing background within 
illuminated regions remains roughly constant. This con- 
trasts with the situation during hydrogen reionization, 
where the same quantity increases in proportion to the 
characteristic bubble size. The difference occurs because 
so few sources contribute to helium reionization: only a 
few are needed to reach the mean post-reionization value 
inside a bubble. During hydrogen reionization, on the 
other hand, the number of visible sources increases pro- 
portionally to the volume of the bubble. 



Despite this smooth match onto the post-reionization 
Universe, there are fairly sharp observational probes of 
the reionization era. Most importantly, we have shown 
that the fraction of pixels with large values of rj (> 1000) 
increases faster than linearly as xhcII increases. Not only 
is the hardness ratio large inside of regions that have not 
yet been fully-ionized, but ionized bubbles without any 
active sources also rapidly recombine and become nearly 
opaque in He II. This should make He II reionization 
easier to identify. 

This contrasts with the post-reionization Universe, 
where rj < 200 everywhere in our model; this maximum 
value corresponds to the minimum ionizing background 
generated by distant quasars. Our model shows that this 
should be a rather sharp cutoff (at least without shadow- 
ing and radiative transfer: iTittlev fc Me iksin 2007), be- 
cause the probability of having zero bright sources within 
two attenuation lengths is very small. Direct compar- 
isons to data are difhcult, however, because real obser- 
vations inevitably contain a number of lower limits to 
r; (corresponding to points where the H I optical depth 
is too small to measure). A careful calibration to more 
detailed simulations is required to draw firm conclusions 
about the timing of helium reionization from the 77 dis- 
tributions. 

On the other hand, the high- J tail of the distribution 
(corresponding to small rj) - which easily reaches val- 
ues several times the mean - remains more or less intact 
throughout at least the latter half of reionization. Thus, 
we expect that measurable pockets of transmission in 
the He II Lya forest will persist into the reio nization 
epoch at z > 3: a true iGuim fc PetersonI (|1965D trough 
will only appear when nearly all of the helium atoms 
are singly ionized. This implies that ultraviolet spectra 
of higher-redshift quasars will continue to provide de- 
tailed information about this important transformational 
epoch, even at redshifts well beyond the nominal "end" 
of reionization. Fortunately, dozens of promising tar- 
gets have now been detected () Zheng et al.ll2004al l2008l: 
ISvphers et al.l [20081 ). and the recent installation of the 
Cosmic Origins Spectrograph and repair of the Space 
Telescope Imaging Spectrograph on the Hubble Space 
Telescope has provided two powerful instruments for such 
exploration. 

Note that this also contrasts with hydrogen reioniza- 
tion, where variations in the i onizing background are 
small during the final phases (jFurlanetto fc Mesingeil 
|2Q09| ). There, the huge number of galaxies responsible 
for reionization markedly decreases the importance of the 
high-J tail. Without such strongly-ionized regions, and 
with the high mean densities at these redshifts, the H I 
Lya forest becomes completely saturated even inside of 
large ionized b ubbles before reion i zation completes; the 
appearance of a lGunn fc PetersonI (jl965l ) trough at z 6 
may or may not hint to us that reionization is ending at 
that time, but in either case the prospects for detect- 
ing substantial transmission at higher redshifts are dim. 
That is not the case for helium, for which our model 
strongly predicts that regions of transmission will be rel- 
atively common throughout the bulk of that process. 

Careful study of the He II forest during the reion- 
ization era may help to illuminate some of the physics 
of that time. For example, comparing the hydrogen 
and helium absorption lines may help us to identify 
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when thermal broadening is important to the forest 
(and hence when, and where, heliuni reionization pho- 
toheats the IGM: iGles er et al.' '200 51: iFurlanetto fc OhI 
I2008at iMcQuinn et all 2009: Bolton_eLalJ [20081). The 
spatial scales over which the hardness ratio fluctuates 
will also constrain the growth of He III bubbles, and of 
the attenuation length. These will no doubt be difficult 
measurements, but they will provide more direct infor- 
mation about the equation of state of the IGM than any 
other method. 

Our model does not, of course, provide a complete pic- 
ture of the ionizing background during helium reioniza- 
tion. Most importantly, we only consider the effects of 
radiative transfer in the crudest possible manner (via a 
fixed attenuation length tq). Shadowing, fluctuations in 
the mean free path, and the large range of relevant fre- 
quencies and spectral indices will all modify the simple 
picture here (IShull et all 12004 iMaselh fc Ferraral 120051 : 
iTittlev fc Meiksinll2007t ). Another aspect tha t we ignore 
is qu asar clustering, which is fairly strong fShe n et al.l 
l2007t iFrancke et al. 2008), altho ugh still probably sub - 
dominant to random fluctuations (jMcQuinn et al.l(2009D . 

Another important detail that our model does not di- 
rectly describe is the evolution of tq. We have argued 
that the attenuation length as seen by a single quasar 
does not evolve significantly (Fig. [4|); however, once a 
universal, more spatially constant background has built 
up, ro will converge to a more uniform value. Reassur- 
ingly, these two appr oaches give similar mean free paths 
at z w 3 (see also IFurlanetto fc OhI l2008bl ). roughly 
where they must converge to each other at the end of 
reionization. On the other hand, in our figures we have 
kept redshift constant for a straightforward comparison; 
if ro evolves rapidly there will be more evolution in /(J) 



than we have shown here (although this is only important 
once the characteristic bubble size exceeds tq). 

In this paper, we have focused on how fluctuations in 
the helium-ionizing background affect the hardness ra- 
tio of forest lines. Of course, the wide variation in J 
even after helium reionization ends may be observable 
in other ways. For example, it will affect the evolution 
of the mean optical depth in the He II forest: strongly 
illuminated regions will allow substantial transmission in 
some reg ions even when most of the Universe is opaque 
(see, e.g. . IFurlanetto fc Mesinge3200S for a similar treat- 
ment for hydrogen reionization). This will cause a more 
gradual evolution in the overall optical depth t^s than 
predicted by naive reionization models. Moreover, the 
high-energy radiation field can also be measured with 
metal lines. The data so far are controversial: some 
measurem ents imply a hardening in the radiation fiel d 
at z - 3 (|SongailalTl998l . l2005l : lAgafonova et al.ll2007f) 
while others are consisten t with no evolution ()Kim et al.l 
l2002t lAguirre et al.ll2004f ) . The complexities of radiative 
transfer during helium re ionization may account for these 
()Madau fc Haardlj I2008D . Alternatively, the strongly 
fluctuating ionizing background throughout reionization 
provides a natural explanation for these differences, al- 
though of course a quantitative estimate is still required. 
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